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Abstract 

Renormalization of composite three-quark operators in dimensional regularization is complicated by the 
mixing of physical and unphysical (evanescent) operators. This mixing must be taken into account in a 
consistent subtraction scheme. In this work we propose a particular scheme that allows one to avoid the 
necessity of additional finite renormalization and is convenient in QCD applications. As an illustration, we 
calculate the two-loop anomalous dimensions of local three-quark operators in this scheme. 



1. Introduction 

The studies of hard exclusive processes includ- 
ing baryons will constitute a significant part of the 
experimental program at the planned new facilities 
at JLAB (Newport News) and FAIR (Darmstadt). 
QCD description of such reactions and also heavy 
hadron weak decays with light baryons in the final 
state (LHCB, CERN) involves baryon wave func- 
tions at small transverse separations, the so-called 
distribution amplitudes (DA) . The nucleon DAs al- 
ready have received considerable attention in the 
literature, see e.g. Refs. [l|, y, [sj, Ij] . In particular 
lattice calculations can provide one with reliable es- 
timates of the lowest moments of baryon DAs which 
are defined by matrix elements of local three-quark 
operators [5|, Isf. 

As usually, the QCD matrix elements depend 
on the scale and on the renormalization scheme. 
The renormalization of baryon operators in dimen- 
sional regularization beyond the leading logarithms 
involves some subtleties that, to our knowledge, 
have not been treated in the literature in a sys- 
tematic way and the purpose of our letter is to fill 
this gap. Our study was fuelled in particular by 
the need to calculate the NLO MS/MOM scheme 
conversion factors in lattice studies (cf. [T], |8|) and 
the NLO extensions of light-cone sum rule calcula- 
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tions of baryon form factors [9| which have already 
started [lO|. 

It is well known that for a generic operator the 
widely used MS prescription does not fix a renor- 
malization scheme completely. Indeed, for a nonin- 
teger dimension d the Lorentz group 5*0(1, d — 1) 
becomes infinite dimensional. As a consequence the 
tensor content of the d— dimensional theory is richer 
than that of a strictly four dimensional one. In 
particular, it means that there exist operators in 
d— dimensions which have no counterparts in four 
dimensions, e.g. totally antisymmetric tensors of 
rank n > 4. Such operators vanish in d = 4 and are 
conventionally referred to as evanescent operators. 
However, one cannot simply exclude evanescent op- 
erators from consideration since under renormaliza- 
tion they mix with the physical operators. The situ- 
ation was thoroughly analyzed by Dugan and Grin- 
They have shown that one can always 
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get rid of the mixing by a suitable finite renormal- 
ization. A more detailed discussion in [12| shows 
that such finite renormalization is not unique. 

We are interested in the renormalization of the 
local three-quark operator 



Qafci 



a/37 






(1) 



Here a, /3, 7 are the spinor indices, while ijk and abc 
are the color and fiavor indices of the quark fields, 
respectively. The flavor structure of the operator 
will be irrelevant for the further discussion, so that 
from now on we will suppress the flavor indices. 
In typical applications one usually tries to get rid 
of multiple spinor indices by contracting ([T]) with 
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suitable 7— matrices. An example is provided by 
the so - called loffe current [13| 

Vj{x) = e'^'= [u\x)Cj^u^{x)] j^rd^ix) , (2) 

or the leading twist nucleon operator 

tjn{x) = e'^'' [u\x)C^u^{x)] ^^{x) . (3) 

Here z is a light-like auxiliary vector and C is the 
charge conjugation matrix. The two-loop anoma- 
lous dimension of the lofFe current was calculated 



14j . Evanescent operators do not contribute at 



this order to the current ^ due to special cance- 
lations however they should be taken into account 
for the current (j31). 

The renormalization of the currents ©, (|31) is 
not of academic interest only |^. The sum rules in- 
volving these currents provide an important tool 
for quantitative study of hadronic properties, see 
e.g. Refs. 0, i, [H, [i3|- The matrix elements of 
the operator ((T|) and the corresponding interpolat- 
ing currents satisfy at tree level certain symmetry 
relations which are derived with help of Fierz trans- 
formations tSj. Therefore it is preferable to use a 
renormalization scheme which respects these rela- 
tions. 

In our opinion the scheme developed in ll|, [l2 , 
14l | is not best suited for renormalization of the 
three-quark operators. We will suggest a different 
scheme and calculate the two-loop anomalous di- 
mensions of the operator ([1]) in this scheme. 

The paper is organized as follows: In Sect. [21 we 
briefly review the renormalization approach taken 
by Dugan and Grinstein [11]. Then we propose a 
new renormalization scheme and discuss its proper- 
ties. In Sect. [3] we present the results of the calcu- 
lation of anomalous dimensions of the operator ([1]) 
in two loops. Conclusions are presented in Sect. SI 

2. Renormalization schemes 



Let us at first discuss the renormalization of 
the loffe current in the scheme [11[. It will be 
more convenient to work with the current f]i [x) = 
IbViix) which definition does not involve the matrix 
75. Next, we will assume that the d— dimensional 
charge conjugation matrix C satisfies the defining 



^The maximal hclicity three-quark operators in the 
AT = 4 SUSY were studied in Refs. [H [3 in AdS/CFT 
context. 



relation C^^C^^ — — 7J. Then one can easily ver- 
ify that the operator fji mixes under renormaliza- 
tion with the operators 



Vn 



Jjk 



u'C-i^'^^u^ 



7 



(n)^k 



(4) 



odd. Here 7^"^ 



„(») 



7 



(") 



7a'" ..A'>^ ® 7m" .-Mn and 



7^1. ../i„ is the antisymmetrized product of n gamma 
matrices, 7^"'..;.„ = "t[^, . . . 7^,.] . 

The renormalized operators fjn take the form|j 



Vln] = ^ ZnkVk 



(5) 



fc=l,3. 



Here and below the square brackets will denote a 
renormalized operator. 

In d— dimensions all operators ^ are indepen- 
dent and satisfy the standard RG equation 

{[MDm + P{as)^o.^5nk + Ink) [Vk] = . (6) 

Here M is the renormalization scale, P{as) is the 
beta function and 7„fc is the anomalous dimension 
matrix 



7 ^-M (-Aziz- 



Z = ZZ-^ (7) 



with Zq being the quark field renormalization con- 
stant, go = Zqq. 

Generic correlation functions with insertion of 
the operator ^„, n > 3, vanish at tree level in d = 4, 
thus these operators are evanescent ones. However, 
it does not hold anymore beyond the tree level ap- 
proximation. Moreover, even if the evanescent op- 
erators were zero at one scale they would reappear 
at another scale due to mixing with physical oper- 
ators. 



It was shown in llj that making a finite renor- 
malization one can ensure vanishing of evanescent 
operators beyond tree level at any scale. In other 
words this means that in an arbitrary scheme the 
evanescent operators 77„, n > 3 can be expressed in 
d = 4 in terms of the first two, "physical" , operators 
fji and ?73 U 



Vln] = ani{as)[i-ii\ +a„3 (««)[%], n > 3. 



(8) 



^In the MS scheme the renormaUzation matrix Z is a 
series in 1/e. 

^ Strictly speaking one can choose any other two oper- 
ators as the basic (physical) ones, however in this case the 
expansion coefficients a„fe will be singular in as- 



The above identity should be understood as an 
equahty between arbitrary correlation functions of 
the operators on the l.h.s and on the r.h.s of this 
equation. (To avoid misunderstanding we remind 
that taking the limit d ^ A one puts 7*^"^ = for 
n > 4 in the renormalized correlation functions.) 

The proof of Eq. ([5|) given in Ref. (ll| employs 
combinatorics of the i?— operation. It was done 
on the example of the four-fermion operators but 
the argumentation is quite general. Another way 
to verify Eq. (jS]) is to note that in some nonmini- 
mal schemes, e.g. the MOM— scheme, the condition 
[r]n\ = for n > 3 holds automatically. Operators 
in any two scheme are related by a finite renormal- 
ization that implies Eq. ([8]). 

Inserting ^ into © one derives the RG equation 
involving the physical operators 771,3 only 

{\^MdM + P{as)d^^5nk + Ink) [flk] = , (9) 

where the indices n, k take the values 1, 3 and 



Ink = Ink + 



/ ^ ^np^pk ■ 
=5,7... 



(10) 



Evidently one has to consider the evanescent oper- 
ators in order to obtain the correct result for the 
anomalous dimension matrix. 

It is also clear that if the operators ??„ enter the 
OPE of certain currents in d— dimensions, the four 
dimensional OPE can be written in terms of the 
physical operators only with the coefficient func- 
tions being appropriately modified. 

In this work we propose a different MS scheme 
for the renormalization of the three-quark opera- 
tor ([T]) (see also the recent paper [13]). Our anal- 
ysis follows closely the approach taken in [19]. In- 
stead of contracting the operator ([ij with different 
7— matrices we will consider the operator with open 
spinor indices. The divergent part of any diagram 
contributing to the correlator 

{Qa(3'y{0)qa'{k)q(3'{p)qY{q)) (11) 

after subtracting the sub-divergences can be cast 
into the form 

/ , Jnmk\£) (1 nmkJajj^ j \^^) 

nrrtk 

where the functions /nmfc(e) are series in 1/e. The 
gamma matrix structures Tnmk sue defined by 

(r„™.)i'!,'^' = 7i">7S^^7S^ (13) 



where it is assumed that all Lorentz indices of 
gamma matrices (jj^...^^) are contracted between 
themselves Q We define the subtraction scheme 
by removing all singular terms (1121) from the cor- 
relator (jlip . Thus the renormalized operator Oap-y 
takes the form 



[OaPj] - Z2p^^ Oa'l3'Y, 



(14) 



where 



■^"/aV " ^+ $I""™'=(^nr«mfc)a/7'' (^^) 



and 



a_,(e)=Ve-Pail,(a.). (16) 



p=i 



The RG equation for the operator [Oap^y] reads 

(^MOm + /3(a,)9„,) [0„^^] = -7:S'''[Oa'/^'7'] : 

(17) 

where the anomalous dimension matrix 7 is given 
by Eq. ([T]). Calculating the inverse matrix Z^^ 
one has to carry out all gamma matrix algebra 
in d— dimensions and this results in emergence of 
e— regular contributions in Z~^. Thus the matrix 
Z~^ contains both singular and regular terms in e 

iZ-'rJ^''' = 1 + E "-'^W i^nmkTj:,''' , (18) 



with a„mk{e) = E^-00 ^ ^ aim/c(«s) • This is dif- 
ferent from the standard situation where Z~^ is a 
series in 1/e, Z^^ — X^fc '^fc(l/^'^) ^^'^ there are no 
finite in e terms. 

The anomalous dimension matrix 7 in 
d— dimensions takes the form similar to (|15|) 






7 — / , ^nmky^s: ^) -L nnik i 



nmk 



(19) 



(we omitted the spinor indices for brevity). The 
coefficients jnmk (cts , e) are regular functions of a^ 
and e. In d = 4 one can drop all F— matrices which 
vanish in four dimensions from the sum (1191). i.e. 



n,m,k < 4. 



^ There is only one nontrivial way to contract all indices. 




Figure 1: The one-loop correction to the operator C'cI/3^ 



Let us see what happens with evanescent oper- 
ators in this scheme. We define the renormahzed 
current [fjn] as follows 



[^„] = (C7("))o;3®7S[£^^''<"K 



(20) 



Since the 7— matrix structure is convoluted with 
the renormalized (finite) operator, one can safely 
put 7^") = for n > 5 in d = 4. Thus, the evanes- 
cent operators automatically vanish in d = 4 in this 
scheme and we can treat the renormalized operator 
[Oa/37] as a pure four dimensional object. Since the 
7— matrices in Eq. (|20p are already four dimensional 
one avoids the problem of defining the matrix 75 in 
d— dimensions. Moreover, this scheme is obviously 
consistent with the Fierz identities, as all contrac- 
tions of the renormalized operator with 7— matrices 
are done in four dimensions. 



The expression for the anomalous dimension ma- 
trix (O takes the form 



2al eZ( 



i)+a(. 



2Z(2)-Z(1)Z(1) 



6oZ(i)) \+0{a''), 



(23) 



where 



Z = ZZ-^ = l + a Z^^l + a^Z^^) + 0{a^). (24) 

Calculating the one-loop diagram shown in Fig. [T] 
one gets the following expression for the renormal- 
ization matrix Z^^) (Z = I + Y,a^Z^''^) 



6e 



r220 + r202 + ro22 + 12 



(25) 



where 



r220 = tE^ ® l^a} ® I 



(26) 



and similarly for the others. The one-loop anoma- 
lous dimension matrix (|23p takes the form 



3. Two loop analysis 



220 



202 



022 



(27) 



In this section we present results of the two-loop 
calculation of the anomalous dimension matrix for 
the three-quark operator ([T]). All calculations were 
done in Feynman gauge. The quark field renormal- 
ization constant Zq in this gauge reads [20| 



In order to present the results of the two-loop calcu- 
lations in a compact form we introduce short-hand 
notations for combinations of 7— matrix structures 
which appear in the calculations 



2^^^ + y 



^(4^^ 



CfCa 



\cJ-^CA^\c, + \Nf 



(21) 



Here a = as/(4^), Cp = {Nj - 1)/27V„ Ca = N, 
and Nf is the number of flavors. The beta function 
for the coupling a is 

Pa{a) = M^ - -2ea + 2boa^ + ©(a^) , (22) 



where 



-2 = r220 + r202 + ro22 , 
'4 = r440 + r404 + ro44 , 

'42 = r422 + r242 + 1^224 : 



r440 -T^f 



^J,l^Jl,2^i3^J■i 



Ta22 =7^''^ 



(4) 



'7Ai 



MlA'2M3/i4 



/, 



'^1/^2 -^ 'M3/i4 



(28) 



where 60 = 2A^//3 - ll/3A^c 



The contributions to the matrix Z from the two- 



loop diagrams read 



in d 



^(2a) 



2nr 
9 U2 



1 



2e 



^(2c) 



^(2d) 



:2 + i2 

1 f 5 / 1 17 \ Nf f I 1 
5 



3 1 4 V £2 30£ ; 6 V £2 6e 



-I- 

36 1 £2 



12 



2e 



3 f 1 



4 [e2 
1 5 



C2 + I2 



.lie.-. 



^""'Tliil'^*-'''^"-^''}- 



^""4{?(5^' + «^= + »» 



1 /I 



£ V4 



I2C2 + 120 



^<->4{^(c. 



1 /I 



£ V2 



(29) 



Here the factor Z(2'^) (Z'^2'')) comes from the quark 
(gluon) self-energy corrections to the one-loop di- 
agram shown in Fig. [TJ Similarly, Z^2c) ^j^(j 2'(2'') 
correspond to the QED - like and QCD - like ver- 
tex corrections, respectively. Z'^2e) ^j^j ^(2/) result 
from to the diagrams with intersecting and non - 
intersecting gluon lines, respectively. Finally, Z'^2ff) 
takes into account contributions from the diagrams 
with all three quarks interacting Q 

The calculations are straightforward so we will 
not present them. Some useful identities for the d- 



11 



dimensional 7-matrices can be found in Refs 
2l| . We also checked that our results are consis- 
tent with the two-loop calculations of Pivovarov 
and Surguladze [ij]. 

Substituting ^ and ^ into ^ and taking 
into account the identity 



C4 



2-2 



:4-2(d-3)C2-3d(rf-l) 



one obtains for the anomalous dimension matrix 



a . 



3 V6 

6a^{bo-l) + 0{a^). 



5C4 - 2{bo - 36)C2 



(30) 



To find the eigenvalues and eigenvectors of the ma- 
trix 7 let us note that the matrices C2 and C4 enter- 
ing (|50| are Lorentz invariant operators. An opera- 
tor O^^^' transforming according to the irreducible 
representation of the Lorentz group which is labeled 
by two spins {j,j) diagonalizes the matrix C2, i.e. 
C2 O(J^) = ci"^ O'-^'^l The corresponding eigenval- 
ues are 



.(1/2,0) _ (3/2,0) _ 



3c 



(1,1/2) _ 



12. 



(31) 



and 
form 



„(w) 



JJj) 



The matrix C4 in d = 4 has the 



C4 == 24(75 75 (g) / -I- 75 (g) / (Ki 75 + / «) 75 «) 75) . 

Its eigenvalues depend on the chirality of the quark 
fields entering an eigenoperator 0^^^\ Namely, 
cj = 3 • 24 if all quark fields have the same chi- 
rality, and C4 — —24 otherwise. 

Thus an eigenoperator carries three quantum 
numbers |j: Lorentz spins (j,j) and "chirality ±". 

In the following we present four operators fea- 
turing different quantum numbers. These are the 
standard operators which appear in the studies of 
the nucleon matrix elements [2|, y, |9|, [isl, |l7| . Two 
of these operators are of twist three and can be 
written in the form |22| 






1 



(32) 



where qL{Fi.) = ^^(1 T 75)9 are left(right)-handed 
spinors and z is a auxiliary light-like vector, z'^ = 0. 
The other two operators have twist four and can be 
written as follows 14 1 



0^^^"^ =£*■'■'= (gfCgi)g^ 



.(^,0) _^»jfe/ .T, 



Ol^-=e*,'-(gf Cgi)4. (33) 



Converting the operators ([5^ . ([55]) into the stan- 

(- 0) (1 -) 

dard notation one finds that OJ ' and 0_ ' ^ cor- 
respond to the lofFe current rjj, Eq. ([2]), and the 



^The diagram where all quark lines connected via triple 
gluon vertex vanishes because of the color structure. 



^The anomalous dimensions are insensitive to the flavor 
structure of the operator so that we disregard it. 



leading twist nucleon current t]N: Eq. ^, respec- Acknowledgments 



tively. The operator O^ 
current [231 



(io) ; 



is related to the Dosch 



■noix) 



— Jjk 



[u\x)C<j^,u^{x)\ 75fT'^"d^x) (34) 



.(§■0) 



whereas 0_,_^ ' belongs to the baryon decuplet cur- 
rent defined in [J]. 

For the anomalous dimensions of the opera- 
tors ((5^ . (|55|) we obtain 






-M 



1 /Of, 

3 

\ TT 
V TT 



2 /9 

4 12 

2 /23 7 

36 " 18 



6o 



,(io) ^ 



TT 
TT 



—A f' 
"12 Y' 



(35) 



The anomalous dimensions of the operator sub- 
set ([33)) were calculated in [14] in another scheme. 
The operators in these two different schemes are 
related to each other by a finite renormalization 



Ia-fO(a^))(( 



0(!'°) = 1 - ^ a -f 0{a^) (oi^^°)) - . (36) 



It can be checked that this factor accounts for the 
mismatch of the anomalous dimensions in the two 
schemes. 

4. Summary 

We proposed a simple scheme for the renormal- 
ization of local three-quark operators in dimen- 
sional regularization. An attractive property of this 
scheme is the guaranteed vanishing of the evanes- 
cent operators in d = 4 dimensions so that one 
can work with physical (four dimensional) opera- 
tors only. The renormalization procedure maintains 
explicitly all identities for operators based on Fierz 
transformations. 

We have calculated the anomalous dimension ma- 
trix at two-loop order and found its eigenvalues. 
Our results for the twist four operators agree with 
the known in literature. The results for the leading 
twist three operators are new. In particular, the 

(1 i) 

anomalous dimension 7_ ' ^ determines the scale 
dependence of the nucleon wave function at ori- 
gin, fN- 
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